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Abstract 

We analyze large phase transitions for U (N) g-deformed two-dimensional Yang- 
Mills theory on the sphere. We determine the phase diagram of the model and we 
show that, for small values of the deformation parameter, the theory exhibits a phase 
transition which is smoothly connected to the Douglas-Kazakov phase transition. For 
large values of the deformation parameter the phase transition is absent. By explicitly 
computing the one-instanton suppression factor in the weakly coupled phase, we also 
show that the transition is triggered by instanton effects. Finally, we present the solu- 
tion of the model in the strongly coupled phase. Our analysis suggests that, on certain 
backgrounds, nonperturbative topological string theory has new phase transitions at 
small radius. From the point of view of gauge theory, it suggests a mechanism to 
smooth out large A^ phase transitions. 
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1 Introduction 

In certain cases, nonperturbative completions of string theory can be obtained by 
considering a holographic description in terms of a D-brane gauge theory. Recently, in 
the case of topological strings, Ooguri, Strominger and Vafa [2Z| made a proposal for 
such a nonperturbative completion based on the connection to the black hole attractor 
mechanism. According to j2Z], the nonperturbative description of topological string 
theory on a Calabi-Yau background is encoded in a D-brane gauge theory living on 
some appropriate cycles of the manifold. 

In |4m 12] this proposal was made more concrete by considering Calabi-Yau back- 
grounds of the form 



where is a Riemann surface of genus g and Li, L2 are hne bundles such that deg(Li) + 
deg(L2) = 2g — 2. In this case, the relevant D-brane gauge theory reduces to a q- 
deformed version of two-dimensional Yang-Mills (YM) theory on the Riemann surface 
Sp. g-deformed 2d YM can be regarded as a one-parameter deformation of the standard 
2d YM theory. As we will explain below, the deformation can be parametrized by a real, 
positive number p, in such a way that as p — 00 one recovers the standard YM theory. 
The g-deformed theory is exactly solvable and one can compute its partition function 
on any Riemann surface. This partition function has a strong coupling expansion as a 
sum over representations of the gauge group, which can be written, following |2I], in 
terms of a product of a chiral and an antichiral sector. The perturbative topological 
string partition function, which was computed in [Hj for this class of geometries, is 
given by a certain limit of this expansion in which the antichiral sector decouples. 
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Once we have a nonperturbative description of the theory, it is natural to ask what 
new phenomena emerge in this description and what their imphcations are for string 
theory. For example, in ^3] the fermionic description of 2d YM on the torus was used 
to study baby universes in string theory. 

2d YM theory on the sphere exhibits an interesting phenomenon: as shown by 
Douglas and Kazakov ^^I, there is a large N, third order phase transition at a critical 
value of the area A = tt^ between a large area phase and a small area phase. From the 
point of view of the small area/weak coupling phase, the phase transition is triggered 
by instantons (20] • From the point of view of the large area/strong coupling phase 
and its string description in terms of branched coverings fH^. E], the transition is 
triggered by the entropy of branch-point singularities jSHl- Due to the Douglas-Kazakov 
transition, the large area expansion of 2d YM theory on the sphere has a finite radius 
of convergence [SB]- 

In this paper we will study the possibility of large phase transitions in g-deformed 
2d YM. Since as the deformation parameter p goes to infinity we recover the usual 
theory, it is natural to expect the transition to occur at large enough values of p. In 
fact, our result show that the transition persists for all p > 2, and we find a critical line 
smoothly connected to the Douglas-Kazakov transition of the standard 2d YM theory. 
We also show that for p < 2, in the regime of strong g- deformation, the phase transition 
does not occur. We also perform a detailed instanton analysis which shows that, as in 
the standard YM case studied in [20], the transition is triggered by instanton effects. 

Most of the analysis of this paper are done in the small area phase. In 2d YM theory 
this phase is described by a Gaussian matrix model. In the g-deformed case, this 
phase is essentially described by the Chern-Simons or Stieltjes-Wigert matrix model 
introduced in and studied in PQEHIEII^- This model, albeit complicated, is exactly 
solvable (in terms of, for example, orthogonal polynomials), and this is the underlying 
reason that we can make exact statements about the location of the critical line and 
the instanton contributions. The large area phase turns out to be more difficult to 
handle. In this paper we present some preliminary results and derive the equations that 
determine the full solution (including an explicit expression for the two-cut resolvent). 
We expect the phase transition of the g-deformed theory to be of third order for p > 2, 
since it is smoothly connected to the transition of Douglas and Kazakov, and indeed 
we give indirect evidence that this is so. 

As in the standard 2d YM, the existence of the phase transition in the g-deformed 
version indicates that the large area expansion has a finite radius of convergence. Ac- 
cording to jHHl 12]; this theory provides a nonperturbative description of topological 

^Connections between Chern-Simons theory and g-deformed 2d Yang-Mills theory have been made, 
from a different perspective, in |2] and 
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string theory on certain Calabi-Yau backgrounds. This suggests that the large area 
expansion breaks down in the full topological string theory, and there is a phase tran- 
sition between a small area phase and a large area phase. From the gauge theory point 
of view, our analysis shows that when the g-deformation is strong enough, the model 
exhibits a single phase. This suggests that g- deformations give a mechanism to smooth 
out large N phase transitions. 

The structure of this paper is as follows: in section 2 we briefly review the Douglas- 
Kazakov transition in 2d YM theory. In section 3 we determine the phase diagram 
of the g-deformed theory and we find a line of critical points parametrized by p, for 
p > 2. In section 4 we adapt the analysis of [201 and study the phase transition of 
the g-deformed theory in terms of instantons in the weakly coupled phase. We find 
an explicit expression for the one-instanton suppression factor which indicates that, 
indeed, the transition is triggered by instanton effects. In section 5 we analyze the 
large area phase, which can be encoded by standard techniques in a two-cut solution 
to an auxiliary matrix model. Finally, in section 6, we discuss the implications of 
our results for topological string theory and outline some problems opened by this 
investigation. 

As this paper was being completed, we became aware of the work j^H], where the 
same problem is studied. After submission, another paper appeared [TT] addressing 
the same issues. 

2 A review of the Douglas-Kazakov transition 

2d YM theory is an exactly solvable model. In particular, the partition function of the 
U{N) theory on the sphere is given by a sum over representations of U{N) (see ^Sj 
and references therein) 



where dim R is the dimension of the representation R, Aisa real and positive parameter 
that can be identified with the area of the sphere, and Ci{R), C2{R) are the first and 
second Casimir of R. We will represent i? by a set of integers {/i, l2, - ■ ■ , In} satisfying 
the inequality 





R 



oo > h > I2 >■■■> In ^ —00. 
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In terms of these integers, the Casimirs have the expression 



N 



(2-3) 



1=1 



Although the above partition function looks rather simple, this theory turns out to 
have a very rich structure. In ^3 1^ it was shown that at large area the partition 
function ()2.H) admits a string representation in terms of branched coverings of Riemann 
surfaces (see [T^ for an excellent review). Douglas and Kazakov found that the planar 
free energy on the sphere exhibits a third order phase transition at the critical value 

A, = tt\ (2.4) 

This large N transition is a continuum analogue of the Gross- Witten-Wadia phase 
transition for 2d YM theory on the lattice [221 E]- Since in this paper we will be 
considering a generalization of the Douglas-Kazakov phase transition, we will briefly 
review how this transition is found. For the rest of this section we will set 6 = 0. 
At large it is natural to introduce a distribution of Young tableaux 

n{x) = j^, x=^- (2.5) 

Defining the shifted distribution 

h{x) = —n{x) + a; — -, (2.6) 
one finds that the planar free energy is given by 

F,{A) = -SgH (2.7) 

where the functional Sclh] reads 

A 3 



Sclh] = — I dx I dy\og\h{x) — h{y)\ + — [ dxh{x] 
Jo Jo 2 Jo 



2 

24 2' 



(2i 



Let us now introduce the density function 



Pih) ^ |. (2,9) 

which is normalized to unity, 

dhp{h) = l. (2.10) 
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One crucial observation of [T^ is that, because of the inequahty ()2.2|) . this density has 
to satisfy 

p{h) < 1 (2.11) 

for all h. We can now write (j2.8p as 

Sg[p] = - J dh J dh'p{h)p{h')log\h~h'\ + ^ j dhp{h)h^-^-^. (2.12) 

This is (up to the p-independent terms) the saddle-point functional for a Gaussian 
matrix model with 't Hooft parameter t = 1/A. It then follows that the density p{h) 
is given by Wigner's semicircle law, 

PG(A,t) = ^v/4t^, (2.13) 

and we find 

p{h)=pG{h,l/A). (2.14) 

However, it is clear that this solution can be valid only for A < tt^, since after this 
point the inequality ()2.11|) is violated. This indicates that there is a phase transition 
at the critical value (|2.4p . 




Figure 1: This figure shows the density p{h) before and after the Douglas-Kazakov 
transition. The solution for A > tt^ can be interpreted as a two-cut solution of an 
auxiliary matrix model. 

For A > tt"^ the Gaussian solution is no longer valid, and Douglas and Kazakov 
argued that one could obtain a solution for the large area phase by considering a a 
density of eigenvalues of the form, 

p(/.) = ^~(^)' -<h<-b,b<h<a, ^2.15) 
1, —b<h<b. 
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where b < a are points in the real positive axis. From the point of view of the density 
p{h), the Douglas-Kazakov transition can be represented as in Fig. ^ ioi A < tt^ the 
Gaussian density gives a good description, but as A > vr^ one finds a new density of the 
form ()2.15|1 . It is easy to see that finding p{h) amounts to finding a two-cut solution 
for a modified matrix model with a logarithmic potential. The explicit solution to 
this problem was worked out in and this allowed them to verify that the phase 



transition at A 



2 ; 



is of third order. It was also verified that the large area solution 



agrees with the string expansion of 

The mechanism behind the Douglas-Kazakov phase transition was further eluci- 
dated in EDI In particular, it was shown by Gross and Matytsin in [20] 
that the Douglas-Kazakov phase transition is driven by instantons. The small area 
phase is dominated by the perturbative vacuum, and instantons are suppressed with 
an exp(— A^) factor. The one-instanton suppression factor at leading order in was 
computed in to be given by 



exp 



-^7gm(^) 



where 



7gm(^) = 27rVvr2 - A - A log 



TT 



+ Vvr2 - A) 



21 



A 



(2.16) 



(2.17) 



Since 7gm(^ = tt^) = 0, as we reach the critical point instantons are not anymore sup- 
pressed and they trigger the phase transition, which is then a consequence of exp(— A^) 
effects which are not visible in the 1/A^ expansion. 



3 The phase diagram of g-deformed 2d YM 

The g-deformed two-dimensional Yang-Mills theory arises as a natural deformation of 
the usual model. This model has been considered in jTOl EHI and more recently, in the 
context of topological string theory, in j2]. The partition function of the g-deformed 
theory on the sphere can be obtained by replacing the dimensions of representations in 
()2.1|) by their quantum counterpart, in the sense of quantum group theory. The result- 
ing partition function depends on the rank N of the gauge group, two real parameters, 
p,gs, and an angle 9. It reads, 

= 5^(dim,i?)'gP^^(^)/2e^^^i(^), (3.1) 

R 

where the quantum dimension of R is given by 

dto,fl= n "•"i^:^r'' . p^) 

l<j<j<Ar J 
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and the g-numbers appearing here are defined as 

[x]=gf-g"i, q = e~^\ (3.3) 
The free energy of the model is defined as 

F' = ^ log Z". (3.4) 
It is convenient to define the parameter A as 

P9s = ^- (3.5) 



As we will see in a moment, A corresponds to the area of the sphere in ()2.1|1 . As in 2d 
YM, we will require A to be positive. Notice that the g-deformed theory is symmetric 
under p, Qg — > —p, —Qs- Therefore, we can restrict ourselves to the range of parameters 
p>0, gs>0. 

An important property of the g-deformed theory is that in a suitable double- scaling 
limit, one recovers ordinary 2D YM. This limit is defined as follows: 

p — > oo, gs — *• 0, A, N fixed. (3.6) 

As — s> with fixed, the quantum dimension becomes the classical dimension: 

dim^i? — >dimR, (3.7) 

and 

(3.8) 

which is the standard weight factor for 2d YM. We then recover the partition function 
()2.1|) for a sphere of area A. The g-deformed theory can then be regarded as a one- 
parameter deformation of 2d YM. 

In this paper we will be interested in the large A^ dynamics of the deformed theory. 
It is useful to introduce the 't Hooft parameter, which is defined as 

t = Ng„ (3.9) 

and we will consider the 't Hooft large A^ limit in which N ^ oo and t and p are fixed. 
The planar free energy 

Fo^(t,p) = lim (3.10) 

N^oo 

will then be a function of t and p. Notice that the limit ()3.6p that gives ordinary 
Yang-Mills theory can be implemented order by order in the 1/A^ expansion by taking 

p — ^ oo, t — ^ 0, pt = A fixed. (3.11) 
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In this way, we recover planar 2d YM on the sphere. We will check many of our results 
for the g-deformed theory, by verifying that in the limit (j3.11|) one recovers the known 
results in 2d YM. 

In order to compute the planar free energy, we follow the steps outlined in the 
previous section for the undeformed theory and represent the planar free energy in 
terms of a functional of a distribution h{x), which is defined as in fl2.(i|l . It is easy to 
see that in the large limit the planar free energy derived from p.ip is given by 

FS{t,p) = -S[h], (3.12) 

where the functional S[h] reads 



S[h] = - dxj^' dy log |2 sinh ^(/i(x) - h{y))\ + ^ 



dxh{x 



,2 



iO / dxh{x) — — + I dx I rfylog |2 sinh -(x — ?/) |, 

Jo 24 Jq Jq 2 



(3.13) 



and in ()3.12j) S[h] is evaluated on the configuration h{x) which minimizes the above 
functional. The last term in ()3.13j] comes from the denominator of the quantum di- 
mension and it is given by 

^ dx rfylog|2sinh ^-{x-y)\ = ^F^'^it), (3.14) 

where 

Fo'^Ht) = ^ - ^ - Li3(e-*) + C(3). (3.15) 



This function is the planar free energy of Chern-Simons theory ^8] , and we recall that 
the polylogarithm of order n is given by 



oo L. 

X 



If we redefine 



the functional ()3.13p becomes 



k=l 



iO 

h{x) h{x) + — , (3.17) 



tp 



S[h]=- f dx f dy log \2 sinh -{h{x) - h{y))\ + — [ dxh{xf 
Jo Jo 2 2 Jq 



pt 0^ 2 , 
l + 2^ + t^^o^(^)- 



(3.18) 
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Since the inclusion of 6 only leads to an additive term in the planar free energy, we 
will set ^ = from now on. After introducing a density function p{h) as in (j2.9p . the 
p-dependent part of the effective action can be written ()3.13|) as 

S[p] = -j dh j dh! p{h)p{h')\og\2sin\i^-{h~h')\ + ^ j dhp{h)h^. (3.19) 

As explained in the previous section, to see if there is a phase transition one first solves 
for the p{h) that extremizes ()3.19|) . assuming a one-cut structure. In order to compute 
p(/i), we have to solve the integral equation derived from ()3.19|1 . 

ph = F j dh'p{h') coth ^(/i - h'), (3.20) 

where P denotes principal value. The density p{h) is supported on a symmetric interval 
(—a, a). A similar integral equation appears in the saddle-point analysis of the Chern- 
Simons matrix model on the three-sphere [33 E In fact, after the change of 

variables j3 = th, ()3.19|1 becomes the planar functional for the Chern-Simons matrix 
model 

^«-/nf n(2^'-">^)^exp{-|f:A0. (3.21) 

i=i i<j ^ i=i 

with 't Hooft parameter ^ = t/p. This connection suggests an effective way of solving 
dSSOl)- As in [UEniEll, we change variables 

A = e^^+'/P, (3.22) 

and we introduce the density for the new variable A, 

P(A) = ^P{h) = ^p{h). (3.23) 
The integral equation ()3.20|) becomes 

^dA'#l. (3.24) 



2t A J A- A 

This is exactly the saddle-point equation for the Chern-Simons/Stieltjes-Wigert matrix 
model, and we can solve it in a variety of ways jH Ell- The direct computation 
performed in jHlI is the most convenient one in view of the two-cut solution that we 
will introduce later, so let us briefly review it. As usual, we introduce a resolvent 

^o(A) = / d\'^^, (3.25) 
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which due to the normahzation ()2.10|) and the redefinition ()3.23|) . satisfies the following 
asymptotic behaviour 

cuo(A) = ^ + 0(A-2), (3.26) 

as A ^ oo. The density p(A) is recovered from the resolvent ljq{X) through the standard 
equation 

pW = -7^ (^o(A + te) - LUoiX - ^e)) . (3.27) 

We are looking for a one-cut solution to the problem, therefore we assume that the 
density of eigenvalues is supported in the interval {a~,a~^), where 

^± = e^t'^+t/p, (3.28) 
It is well known that ujq{X) can be computed as 

where C is a contour around the cut (a~, a"*"), and 

^7(A) = |^, r(A) = v/(A-a-)(A-a+). (3.30) 




Figure 2: This figure shows the deformation of the contour needed to compute the 
resolvent in ()3.29|1 . We pick a residue ai z = p, and we have to encircle the singularity 
at the origin as well as the branch cut of the logarithm, which on the left hand side is 
represented by the dashed lines. 

The standard way to compute an integral like (j3.29|) is to deform the contour. Since 
the logarithm has a branch cut, we cannot push the contour to infinity. Instead, we 
deform the contour as indicated in Fig. |21 We pick a pole at z = A, and then we 
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surround the cut of the logarithm along the negative real axis and the singularity at 
z = with a small circle of radius e. A similar situation appears in, for example, 
The final formula for the resolvent is 



plogA p 



dz 



dz log z 



z{z — X)r{z) Jc^ 2ni z{z — X)r{z) 



.. (3.31) 



The integrals in the second line have loge singularities as e — 0, but they cancel each 
other, and after some computations one finds for the resolvent: 



'^o(A)=-|^log 



+ 



V 



2tA ycra^ 



r(A) 



log 



4a 



(3.32) 



In order to satisfy the asymptotics ()3.2(jj) the second term must vanish, and the first 
one must go like 1/A. This implies 



4a =2vaa+ + a + a^ 



(3.33) 



and from here we obtain the positions of the endpoints of the cut a , a"*" as a function 
of t/p: 

a- =2e2*/P - e*/P - 2et Ve*/p - 1, 
a+ =2e^'/P - e''P + 2ei \fe^J^^. 
The final expression for the resolvent is then 



^o(A) = -^lo^ 



1 + e-*/f A + ^{l + e-*/f A)2 - 4A 



2A 



and from here we easily find the density of eigenvalues 



p 

P(A) = — tan-i 



v/4A - (1 + e-*IP\y 
1 + e-*/PA 



We can now go back to the original variable /i, to find 



p{h) = — tan"""^ 

TT 



e^/p' -cosh\Ah/{2p)) 
cosh{Ah/{2p)) 



(3.34) 



(3.35) 



(3.36) 



(3.37) 
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which has its support on (—a, a) with 

a = ^cosh-i(e^/(2p')). (3.38) 

As a test of this result, notice that in the double-scahng hmit ()3.11|1 one finds 

p{h)=pG{hMA) + 0{l/p^), (3.39) 
therefore the leading term is exactly the Wigner semi-circle distribution obtained by 

mi- 

In order to assess the possibility of phase transitions, we have to verify the condition 
pTT]) . Notice first that |tan-i(x)| < f, therefore 

p{h) < p/2 (3.40) 

for all h. A first conclusion is that there is no phase transition for p < 2. For p > 2 
there is indeed a phase transition which occurs when the value of A is such that the 
maximum of the distribution reaches the value 1. Since the maximum occurs at /i = 0, 
we immediately find the following line of critical points: 

A,{p) = pHog(^l + tan" (^^^y p>2. (3.41) 

As J9 oo, 

A,{p)^TT^, (3.42) 

in agreement with the result of Douglas and Kazakov ()2.4|) . Notice that A*(p) is a 
decreasing function oi p for p > 2, and as p ^ 2^, the critical area increases to infinity. 
For a given p, the small area phase occurs for A < v4*(p), and in this phase the planar 
free energy is well described by the distribution ()3.37p . 

We then have the phase diagram represented in Fig. El The horizontal axis represents 
the parameter p, while the vertical axis represents A. The critical line, described by 
the function 1)3.411) . has two asymptotes, represented by dashed lines: as p — > cxd it 
approaches the horizontal dashed line A = n"^, which corresponds to the Douglas- 
Kazakov phase transition. As p — 2+ it approaches the vertical asymptote. For 
p G (0,2] there is no phase transition. Notice that, if we parametrize the planar 
g-deformed theory in terms of p and A, the region p (yo corresponds to a small 
deformation, while the region p < 2 corresponds to a large deformation. We then see 
that, if we start with ordinary 2d YM and we turn on the deformation parameter 1/p, 
the Douglas-Kazakov phase transition persists although the critical area increases. At 
p = 2 there is a "barrier" where the critical area becomes infinite. Therefore, when the 
deformation parameter is large enough, the large phase transition is smoothed out. 
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Figure 3: This figure represents the phase diagram of g-deformed 2d YM theory. The 
horizontal axis represents the parameter p, while the vertical axis represents A. The 
curve shown in the figure is the critical line (j3.41|) . which separates the phases of small 
and large area. The horizontal dashed line, which is the asymptote of the curve as 
p —>■ oo, represents the A = n'^ critical point of Douglas and Kazakov. 

To find the free energy in the small area phase, we have to compute the functional 
(I3.19|l evaluated on the density ()3.37|1 . Since this functional is closely related to the 
functional describing the planar Chern-Simons matrix model, we can borrow the results 
from ^HllHl]. From [S^ it follows that, at large N, the matrix integral ()3.21|) is given 
by 

exp(iV^Fo(0), (3.43) 

with 

MO = ^2^0^(0 + 4 (3.44) 
Since ^ = t/p in our example, we finally obtain 

m,P) = Y2{p'Po"('/P) - ^Po'it)) + + (3.45) 
As a further check of this expression, notice that, after using the expansion, 

Li3(e-*) = C(3) - + - ^ log + (3.46) 
one finds in the double- scaling limit ()3.11|) 

i^o'(t,P)-^-^logA+^, (3.47) 
which indeed is the free energy of the usual 2d YM theory in the small area phase. 
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4 Instanton analysis 



Since g-deformed 2d YM theory is a one-parameter deformation of the standard one, 
we expect the phase transition discovered in the previous section to be triggered by 
instantons as well. In this section we will verify this by computing the one-instanton 
suppression factor in the g-deformed case. This will also give an intuitive explanation 
of why the phase transition is absent for p < 2. 

The starting point of the discussion is to write the partition function of the theory 
in a way that makes manifest the instanton content of the model. Since g-deformed 2d 
YM theory has the same action as standard 2d YM, but differs in the measure , we 
expect the partition function to be expressed in terms of a sum over instantons. 



Z^ = 5^Mn,)exp(-^f;n^^], (4.1) 



where nj, j = 1, ■ ■ ■ , N, are the instanton numbers characterizing a classical solution 
|2Uj . and w{ni) is the weight of such a configuration in the semiclassical expansion. In 
order to compute the weights w{nj), we follow the technique used by Minahan and 
Polychronakos |^ in standard 2d YM and perform a Poisson resummation of the 
original expression p.ip . This can be regarded as a duality transformation which takes 
us from the large A phase where the expansion (jH.lj) is valid, to the small area phase 
where the semiclassical expansion ()4.ip is valid. The partition function can then be 
written as 

Z'' = Cj]F2(27rn,), (4.2) 
where F^i Fourier transform with respect to the variables pj = Ij — j + 1/2: 



F2{x,)= /]^dp^.e-^E.-.^'.]^('2sinh^(p,-p,)) exp(-A^p2^, 



(4.3) 



and we are setting the Q angle to zero. This transform can be performed by first 
computing 

^i(^.)= / n^P^-^'^^^"^'^n(2si^hl(p,-p,)]exp(-A^p2^^ (4.4) 

and then doing a convolution. The integral ()4.4|) reduces to a Gaussian after using 
Weyl's denominator formula for a general Lie algebra, 

J2 e(w)e"^^^-" = n 2 sinh (4.5) 
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where a are the positive roots, w G W are the elements of the Weyl group, and e(w) 
is the parity of w. We find, up to a multiphcative constant. 



N ^ 



Fi{xj) =exp(^-— Jj2sin— (xj -Xfc), (4.6) 



3=1 j<k 

and using convolution we finally obtain 



F2{xj) =exp(-— 5^ 

X / n n ( ^ ^^"^ ^ (^^'^ + ^jfc) ^^"^ ^ (^^'^ " ^i^) ) (~ ^ 
1=1 ,<fc ^ / 1=1 



(4.7) 

where we introduced the notation Xjk = Xj — Xk- The instanton weight has then the 
expression 



f ^ f t t \ N ^ 

«^(%) = / n ^^J' n ( ^ ^(27™jfc + Vjk) sin ^(2vrnjfc - ) exp ^ yj^' 
i=i j<k ^ ^ j=i 

(4.8) 

which is a g-deformed version of the result in PHI for standard 2d YM. 

As it was pointed out in [20], a precise way to evaluate the importance of instan- 
ton contributions to the partition function is to compare the contribution of the one- 
instanton term in the semiclassical expansion ()4.1|) to the contribution of the pertur- 
bative vacuum. The relative weight of these contributions defines a function '~f{A,p) as 
follows 



exp 



-^l{A,p) =exp(-iV^)— , (4.9) 
A J V A J Wo 



where the exponent in the right hand side involves the instanton action for rii = 
l,'^i>i = 0, and we have denoted 

w, ^ w{10,...,0)_ ^^^^^^ 
Wo w(0, ■ ■ ■ ,0) 

We call the function in (|4.9p the one-instanton suppression factor. Notice that, as long 
as 7(A,p) is different from zero, instantons will be suppressed in the large limit. 
The suppression is bigger the larger 7(A,p) is. In the remaining of this section, we will 
compute '-f{A,p) in the small area phase of g-deformed 2d YM, and we will study its 
properties. 



16 



Let us first define the partition function 



Zn= lldy,ll( 2 sin— iy,~y,)\ exp(-— J^t/I)- (4.11) 

This is very close to the partition function of the Chern-Simons matrix model, although 
it has a sin interaction between eigenvalues instead of a sinh interaction. We can then 
use the results of the previous section after changing 

P^-^f, A^j, (4.12) 

and doing carefully the analytic continuation of p to the imaginary axis. Equivalently, 
we can change variables y = —iA(] /t in ()4.11|) to obtain the matrix model ()3.2H) with 
^ = —A/p^. One can then see from the formulae presented in the last section that the 
planar limit of ()4.11|) is controlled by the following density of eigenvalues, 



Civ) = — 7 tanh""^ 
^^^^ ttA 



V^cos2(?//(2p)) -e-^/p' 
cos{y/{2p)) 



(4.13) 



with endpoints located at 

r = 2pcos-^(e-^/('*'')). (4.14) 

As p ^ CX3, one can easily check that ({y) pciv, A). 

We can now evaluate ()4.10|1 . Notice first that wq = Zn- On the other hand, as in 
PU] . one has 

Wi = 

/oo ( n (4 sin ^ (2vr + (y, - y^)) sin 1 (27r - (y, - y,))^ )^_^, 

(4.15) 

where the correlator is computed in the model ()4.1H) with — 1 variables. Since we 
are interested in the large behavior of the one-instanton suppression factor, we can 
compute the different integrals in the saddle-point approximation. This in particular 
means that we can set yi = inside the correlator in ()4.15|) . We find, 

Wi _ 
Wq 

(^) ^''''^{^^"^^ / ^^C(y)log(4sin^(27r + y)sin^(27r-y))|. 

(4.16) 
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We have first to evaluate the quotient Z^v-i/^at in the large N limit. It is easy to see 
that, at leading order in N, this quotient is 

exp{-iV(2Fo(0+ei^o(0)}, (4.17) 

where Fq{^) is given in (j3.44j) . Here, ^ = —A/p"^, and after an analytic continuation 
^ — s> — ^ we find, 

2Fo(0 + ei^o(0 = ^(Li2(e-^/^^) - ^), (4.18) 

up to an overall sign (—1)^ in Z^^i/Z^. Putting everything together, we obtain the 
following formula for the function 'j{A,p) defined in ()4.9|1 : 

f .1 1 X (4.19) 

-A J dy C{y)\ogi^A sin— {27V + y) sin— {27T -y)y 

The integral in ()4.19j) can be evaluated analytically. Notice first that in any matrix 
model one has 



F{v) = I dAp(A)log(l - X/v) = I dv'{uJo{v') - 1/v'). 

J J oo 



(4.20) 



This follows directly from the definition of the resolvent in ()3.25|) . Taking into account 
the redefinition (j4.12j) . we find that the integral in (j4.19|) is given by 

2ReF(e-^/^''+2"'/P), (4.21) 

where F{v) is obtained as in ()4.20|) . and the relevant resolvent is ()3.35|) . After some 
work, and using standard identities for the dilogarithm, one finds the following expres- 
sion: 

7(A,p) =2vr2 -p2(Li2(e-^/^') + + 2p'ReG{f^{p, A)J.{p, A)), (4.22) 

where ^ 

y) =2*^^°S + log xlog (1 - y) + Li2(l -x)+ Li2(?/), 



f±{p,A) =e^^[±A/{2p')+i{^-7i/p)y (4,23) 
if = tan" 



^ f a/c^^/p^ — cos2(7r/p) 



cos(7r/p) 

In order to understand the properties of the instanton suppression factor, we have 
studied (analytically and numerically) the properties of ()4.22|) as a function of A and 
p for p > 2, A < A^,{p). The main results of this analysis are the following: 
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Figure 4: This figure shows the function 'y{A,p) appearing in the one-instanton sup- 
pression factor, plotted as a function of A, and for the values p = 2.1, 3, oo, from top to 
bottom. For each p it is a decreasing function of the area and vanishes at the critical 
value A^{p). 



1. As p ^ oo, the function 'y{A,p) becomes the function 7gm(^) introduced in 
(I2.17|) . This is a consistency check of the solution. 

2. For any fixed p > 2, 'j{A,p) takes the value 27i^ at A = and then it decreases 
monotonically as the area is increased. At the critical area (j3.4ip one has 

7(A,(p),p) = 0. (4.24) 

The vanishing of 'y{A,p) at the critical area can be proved analytically, since a.t A = 
A(p), 

/±(p, AM) = {cos{n/p)y\-'^-/P. (4.25) 

For arguments of this form (which are algebraic numbers) the dilogarithm satisfies 
nontrivial identities |^ that can be easily shown to lead to ()4.24|) . 

3. For p < p', one has that '~f{A,p) > 'y{A,p') in their common range A < A^{p'). 

These properties are illustrated in Fig. HI which shows the function '~f{A,p) as a func- 
tion of the area for the values p = 2.1, 3, oo, from top to bottom. The above properties 
show that the one-instanton suppression factor in the small area phase decreases as the 
area grows, until it vanishes at A^,{p). Therefore, at the line of critical points found in 
section 3, the instantons are not suppressed anymore and they become favorable con- 
figurations. This shows that the phase transition for the g-deformed theory is indeed 
triggered by instantons, and follows a mechanism similar to the one studied in |2()j : 
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for A > A^,{p), the entropy of the instantons dominates over their Boltzmann weight. 
The above analysis also shows that, as p decreases, the instanton suppression factor 
becomes bigger and bigger, pushing the critical value of the area to ever larger values. 
This indicates that the smoothing out of the phase transition for p < 2 is due to the 
fact that the instantons are suppressed for all values of A and we only have one phase 
dominated by the perturbative vacuum rii = 0. 

5 The two-cut solution 

In this section we give some preliminary results about the large area phase of the 
theory. After the phase transition found in section 3, we expect a distribution p{h) a 
la Douglas-Kazakov, with the shape shown in the r.h.s. of Fig. ^ and characterized by 
two points d, b. The distribution governing the large area distribution is then of the 
form 

rp(/.), -a<h<-brb<h<a, 

^ [1, ~b<h<b. ^ ^ 

After changing variables A = exp{th + t/p) as in the previous section, the new density 
of eigenvalues p(A) = l/{t\)p{h) has support on the two intervals (a~,6~), (6"*",a+), 
where 

a± = e*/P±*^ b^ = e'/P^'K (5.2) 
This density satisfies the following integral equation, 

P logA 1, A/6+-1 ^ /■ P(A') ,rn^ 





b- b. 



Figure 5: This figure shows the deformation of the contour needed to compute the 
resolvent in the two-cut solution. We have to encircle the singularity at the origin, and 
the two branch cuts denoted by thick lines on the left. 
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As in the one-cut case, we introduce a resolvent 



'^o(A) = J YZ^dX'. (5.4) 
This can be again computed by the contour integral ()3.29|) . but now 



r{z) = ^/{z - a-){z - a+){z - b~){z - b+). (5.5) 

and C is the union of the contours surrounding the cuts (a^, b^), (a+, 6+). To perform 
the integral ()3.29p we deform the contours in the way shown in Fig. El we now encircle 
the branch cut along (— cxd, 0), coming from log A, and the branch cut of the integrand 
along {b~,b~^). The answer for the resolvent is 

. plogA 1 A/&+ - 1 

"°^^)=2^^ + a^°^A7F3T 

- '^'^2tH[-j_^ ziz-X)riz) + 1 2^^ziz-X)riz)] ^'''^ 
^r(A) f''^ dz 



t Jb~ z{z-X)r{z)' 

The above integrals can be expressed in terms of elliptic functions. We will change 
notation a~ ,b~ ,b'^, to d, c, b, a. Define 



IiX,u) 



dz 



{z - X)r{z) 

2 r ^ 

[c - rf)n(0, n, k) + {d- A)F(0, k) 



(A - c)(A - d)y/{a-c){b-d) 

where 11(0, n, k) and -F(0, k) are incomplete elliptic integrals of the third and the first 
kind, respectively, and 

sin2 = ^ _ {a~d){X-c) ^ {b - c){a - d) 

{a — d){c — u)' {a — c){X — dy {a — c){b — d) 

In what follows it will be convenient to introduce the following angles 0i, 02 and 
variables rii and 

sm 01 = -, sm 02 = -, 

a — d ca — d 

ni = , n2 



a — c d a — c 
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In terms of these variables one finds, 

/(A, -oo) =- f fr=^S - ^)n(0i, ^, k)^{d- A)F(0i, k) 

(A — c)\ \ — d)\J [a — c)[b — a) 

0) =- J ^ == \ (c - d)n(02, n, k) + id- A)F(02, /c) 

(A — c)[X — d)^^ [a — c)[b — d) 

/(0,-oo) = =J===|(c-rf)n(0i,n2,A;) + c/F(0i,A;) 



cd\J (a — c)(6 — (i) 

(5.10) 

The first integral in the second line of ()5.6|) is given by 

J (/(A, -oo) - /(A, 0) - J(0, -oo) + /(O, -e)) . (5.11) 

The second integral in the second line is simply a residue and it can be computed 
immediately: 

-^^l=loge. (5.12) 
^ yabcd 

We now compute /(O, — e) at next-to-leading order in e. This will have a logarithmic 
singularity which will cancel ()5.12j) . In order to do that, we need the following identity 

n(0, n, k) = 5(n) |- log ^ - ^^1' (5-13) 

where 

^^''^=\{l-nW-n)) ' ^=2^' ^ = 2^(^^ ' ^'-^'^ 

and the r parameter in the theta functions is given as usual by 

q = e^^'^ = exp{-7r K'{k)/K{k)). (5.15) 

Notice that, when 

sin^ = - (5.16) 

n 

we have a logarithmic singularity in the elliptic integral 11(0, n, k). This is immediately 
checked in the integral representation of the elliptic function^. We can now use ()5.13|) 
to extract the next-to-leading behavior. Since 

sm (p = , n = n2, (5-17) 

a — d c + e 



■^Notice that in the conventions we are using the n in 11(0, n, k) corresponds to —n in the definition 
given in |H|. 
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the leading behaviour of n(0, /c, 712) is given by 



'^"^^og. + 0(.), 



(5.18) 



2 

where /32 is given by 



This leads to the following expression 



/(O, -e) = -^i= log e + /(O, 0) + 0(e), (5.20) 
V aoca 



where 

Ca/ (a — c)(6 — c?) ' ^ 
From the above result we sec that the singularities as e — >■ cancel, as wished. 
We now consider the remaining integral. Define 

dz 



(5.21) 



J(A) ^ / 

J c 

^ a-b)U{m,k) + {b- X)K{k) 



{z — X)r{z) 

2 



(A - a)(A - b)^y{a-c){b-d) 

where 



(5.22) 



(b-c)(X-a) 
We then have, 

/(73|;;w4(-'W-^(°))' (^•^*) 

where J(0) is given explicitly as 

^ |(a-6)n(m(0),A;) + 6X(A;)|. (5.25) 



J(0) 



061,/ (a — c)(6 — d) 

Putting everything together, we find the following expression for the resolvent: 
. plogA 1 \/b-l 

pr(A) 



+ ^ (/(A, -00) - /(A, 0) - /(O, -00) + /(0, 0)) (5.26) 

.rgl(.(A)-.(o)). 
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As A — > oo, this is indeed a Laurent series in A: using again ()5.13|) it is easy to see 
that /(A, — oo) contains a term of the form — log(A)/r(A) that cancels against the first 
term in ()5.26p . In order to derive the conditions for the endpoints of the cut, we must 
impose the asymptotic behaviour 

CuoiX) = ^^ + OiX~'). (5.27) 

We find three conditions. First of all, notice that there is a term of order A coming 
from the integrals /(O, 0), /(O, — oo), and J(0). Imposing the cancellation of this term, 
one obtains the condition 

p(/(0, 0) - /(O, -oo)) - 2 J(0) = 0. (5.28) 
The vanishing of the constant term leads to the condition 

p(f(02, k) - F(0i, k)^ = 2K{k), (5.29) 
Finally, the fact that the 1/A term has the coefficient 1 — 2b leads to a third condition. 



pi^{a + b + d - c){F{4>i, k) - F{4>2, k)) - 2{c - d)U{4>2,ni, k) 

. v/raF^(-2ft|||| + >o. ^ - .og(^.(nO) - log 

+ 2{b + d + c- a)K{k) + 2{d - b)Il{m^, k) = t, 

(5.30) 

where 

F(02,fc) b-c 

= -mW' = ^'-'^^ 

These conditions determine the endpoints d,b as functions of the parameters t,p. We 
seem to have three conditions for two unknowns, but since we started with a symmetric 
problem and we just changed variables, one of the conditions is redundant. This is 
not easy to verify from the above expressions, but can be checked, for example, by 
doing a small t expansion of the equations, and assuming a power series ansatz for the 
endpoints: 

oo oo 

d{t,A) = Y,^n{A)r, b{t,A) = Y,bn{A)r. (5.32) 

n=0 n=0 

The ansatz is justified by the fact that, as t — > with A fixed, we must recover the 
standard YM result obtained in ^Jj. One can see that, at leading order in t, the three 
conditions above lead to the same equation, namely 

^^A = 2Kiko), (5.33) 
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where 



kl = (5.34) 
(ao + 60)^ 

Using standard properties of elliptic functions, one can easily check that this condition 
becomes 

A = -K{bo/ao), (5.35) 
ao 

which is precisely one of the equations found in [17\ Notice that, by making use of 
()5.28|) . we can simplify the expression of the resolvent to 

. plogA 1 A/6-1 

— -S-V^ (5.36) 

In principle, the above conditions for a, 6, together with the explicit expression for the 
resolvent in ()5.36|) . determine completely the solution for the large area phase. These 
conditions are rather intricate to be treated analytically, but one could study them 
numerically. 

The most important question to address is the order of the phase transition for 
different values of p. This of course can be seen, as in [T7j, by computing the free 
energy in the large area phase that we have just analyzed. Since the line of critical 
points is smoothly connected to the Douglas-Kazakov transition, we should expect the 
transition in the g-deformed theory to be of third order for any p > 2. Indeed, one 
can find indirect evidence that this is the case by using an argument in [213] based on 
double-scaling limits. If we consider a theory with a large n-th order phase transition 
at a critical area A = A^, between phases I and II, the free energy has the following 
behaviour 

F',{A) - Fo"(A) ~ (A, - Ar. (5.37) 

To define a double-scaling limit of such a theory, one should introduce a string coupling 
constant fig through 

1^;^ = N\A,- A)"". (5.38) 

The nonperturbative effects of such a theory are expected to be of the form exp(— 1 / fig). 
But this means that the instanton effects in the original theory should have the be- 
haviour exp{—N'y{A)), with 

7(A) ~ {A, - Ar/\ (5.39) 

Indeed, in it is found that the function ()2.17|) appearing in the instanton suppres- 
sion factor has exactly the behaviour ()5.39j) with n = 3 near the Douglas-Kazakov 
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transition point, as required for the existence of a double-scaling limit at a third order 
phase transition. According to this argument, the behaviour of the instanton suppres- 
sion factor near the critical point can be indeed regarded as an indirect way to probe the 
order of the phase transition. We have checked numerically that the function 'y{A,p) 
that we found in ()4.22|) behaves indeed as 

7(Ap)~(A(p)-^)'/' (5.40) 

near for various values of p > 2. This is indeed consistent with the large 

phase transition of the g-deformed theory being of third order for all p > 2. 

We should also mention that, in j4j, general criteria have been formulated to de- 
termine the order of a phase transition for a model based on a distribution of Young 
tableaux. These criteria only depend on the behaviour of the density ()3.37|) in the small 
area phase. It can be easily seen that according to these criteria, the phase transition 
of the g-deformed theory is of third order for any p > 2. 

6 Implications for topological string theory and open 
problems 

In this paper we have shown that g-deformed 2d YM theory exhibits an interesting 
phase structure, with a Douglas-Kazakov phase transition smoothly connected to that 
of the standard YM theory, and a "barrier" at p = 2. One of the original motivations 
of this analysis was the appearance of the g-deformed theory as a nonperturbative 
completion of topological string theory on certain Calabi-Yau backgrounds, g-deformed 
2d YM on the sphere has been proposed in |^ |2] as a nonperturbative, holographic 
description of topological strings on the local Calabi-Yau manifold 

Oi-p) ® Oip - 2) ^ P\ (6.1) 

where the integer number p > corresponds to the parameter p appearing in ()3.1|1 . 
Explicit computations in show that the perturbative partition function computed in 
[Oj appears as a certain decoupling limit of the large area expansion of the g-deformed 
theory. However, the fact that this theory exhibits a phase transition suggests that, 
for geometries of the form 1)6.11) with p > 2, the large area expansion has a finite 
radius of convergence which, in terms of the 't Hooft parameter t, is given by t^:{p) = 
^*{p)/P- As p becomes larger, the radius of convergence becomes smaller. Therefore, 
the conjecture of |SZI suggests that for the geometries ()6.1|1 with p > 2, there will be 
a phase transition at small radius in the full, nonperturbatively completed topological 
string theory. What are the possible interpretations of this phase transition in the 
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topological string theory context? We will mention here three possibilities, although a 
better understanding of the implications of the phase transition of g-deformed YM to 
nonperturbative topological strings will require a more detailed treatment^: 

1. A first possibility is that the phase transition in the g-deformed theory indicates 
a topology change in the Calabi-Yau background. After all, the small and the large 
area phases are described by different master fields of the two-dimensional theory, 
corresponding to the one-cut and two-cut solutions discussed above, and it is known 
that in large dualities the master field encodes the geometry of the target PHIC]- This 
topology change might be also interpreted, as in [T3] , in terms of a process involving a 
splitting of baby universes. 

2. A second possibility is that the small area phase does not have a geometric 
interpretation. One indication of that is the string description of standard 2d YM: the 
analysis of 1^ shows that the large area expansion has an interpretation in terms 
of branched coverings of the sphere. However, it has been argued that the existence 
of a large A^ phase transition suggests that this geometric picture does not hold for 
the small area phase |20ij. In the same vein, it is likely that the small area phase 
of the g-deformed theory is not described appropriately by topological strings with a 
geometric target. This is in fact very reminiscent of the analysis of [H] (see also |Z1E1); 
where it was shown that the large A^ phase transition of the unitary matrix model 
corresponds, in AdS/CFT at finite temperature, to the point where the horizon of 
the small AdS black hole becomes comparable to the string scale. At this point, the 
supergravity/geometry picture breaks down. The situation we are considering here 
could be a topological string analogue of the large A^ transition of [H]. 

3. A more conservative possibility is that the conjecture of j2ZI does not fully apply 
to the local geometries ()6.1|) when p > 2, or at least does not apply to the small area 
phase. The original conjecture was formulated for compact Calabi-Yau threefolds, and 
there may be subtleties when applying it to the noncompact case. It turns out that 
precisely for p > 2 there are obstructions for contracting the inside (Kj.lj) to a point 
[HH 5 and because of this reason one can expect these geometries not to arise as a 
decompactification limit of a compact Calabi-Yau. It is intriguing that the "barrier" 
p = 2 that we found in this paper is the same that occurs in the geometric setting. 

In extracting the consequences of our analysis for the nonperturbative physics of 
topological strings, there is another point that should be mentioned. In our analysis 
we considered the saddle-point solution of the functional S[h], and we found that this 
leads to a distribution where (h) = and the dependence on the 6 angle is trivial. 
However, it has been argued in by studying the instanton weight factors, that the 

"^The discussion that follows benefited greatly from conversations with M. Aganagic, D. Morrison, 
H. Ooguri and N. Saulina. 
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presence of a nonzero 6 changes the location of the critical line. This is an interesting 
possibility and deserves further study. Also, we have restricted ourselves to solutions 
with zero U{1) charge. This is indeed the true vacuum of the theory [3], but one 
could also consider saddle-point solutions like those in : one imposes the constraint 
(h) = Q, where Q is the U{1) charge, solves for the density, and then finally sums over 
all integer charges with a weight exp{iQ6). It may happen that, in order to compare our 
results with those of 0, one should use this prescription to include the U{1) charges. 

It is also worth pointing out that the instanton weight factors considered in section 
4 are closely related to the degeneracies of BPS states analyzed in j2]. It is likely that 
the techniques of jSHI that we used and extended to the g-deformed case in order to 
compute these weights lead to a useful technique to obtain the degeneracies. 

From the point of view of the two-dimensional gauge theory, the results of this 
paper indicate that, when the deformation parameter is sufficiently large, the large N 
phase transition is smoothed out already at the planar level. This is an interesting, 
new mechanism for smoothing out large N transitions which may have implications in 
other contexts (the other mechanism we are aware of to smooth out these transitions 
requires performing a double- scaling limit, as in p^lH]. and involves a resummation of 
the 1/N expansion). 

There are also various open questions concerning the gauge theory aspects of our 
analysis. Of course, the two-cut solution that we presented in this paper should be 
investigated in more detail. One could also investigate the phase structure and free 
energy of the chiral version of the g-deformed theory (in the 2d YM case, this has 
been done in ^1 ED!)- Since the chiral sector makes a more direct contact with the 
perturbative topological string amplitudes, this may help in understanding better the 
holographic description proposed in [101 12] . It would be also very interesting to analyze 
the subleading corrections to the planar result in the small area phase. In [20] 
this was done for the standard YM case by using a discretized version of orthogonal 
polynomials, but it is not obvious how to generalize this to a discrete model with a sinh 
interaction. Such a generalization would make it also possible to define a double-scaled 
theory near the critical line of the g-deformed theory, as we briefiy discussed in the last 
section. 
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